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Abstract 

The interpretation of the J/ip suppression patterns observed in nuclear collisions, 
at CERN and RHIC, cLS cL SI gnature of the formation of a deconfmed phase of 
QCD matter, requires knowing which fractions of the measured J/ip yields, in pp 
collisions, are due to decays of heavier charmonium states. From a detailed anal- 
ysis of the available mid-rapidity charmonium hadro-production cross sections, 
or their ratios, we determine that the J/ip feed-down contributions from ip' and 
Xc decays are, respectively, (8.1 ± 0.3)% and (25 ± 5)%. These proton-proton 
values are derived from global averages of the proton-nucleus measurements, as- 
suming that the charmonium states are exponentially absorbed with the length 
of matter they traverse in the nuclear targets. 
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1 Introduction and motivation 



In the very hot and dense, strongly interacting matter produced in high-energy nu- 
clear collisions, it is expected that the QCD binding potential is screened, the screen- 
ing level increasing with the energy density of the created system [1] . Depending on 
the screening level, it may happen that the charmonium states "dissolve" into open 
charm mesons [2]. Since different quarkonium states have different binding energies, 
they are expected to dissolve at successive "thresholds" in the energy density or tem- 
perature of the medium [3]. In particular, the if)' and Xc states should be easier to 
"melt" than the more strongly bound J/if> state. Therefore, a "spectral analysis" of 
the charmonium production yields, in several collision systems (from light to heavy 
nuclei) and in several collision centralities (from peripheral to central), should provide 
very interesting information concerning the nature of the produced matter. 

Experimentally, it has not yet been possible to directly measure the production 
yields of the \ c state in heavy-ion collisions. However, it is well known that a signifi- 
cant fraction of the J/ip mesons observed in pp collisions are, in fact, produced by Xc 
radiative decays. The J/ip production yield measured in heavy-ion collisions could 
thus show a significant level of suppression, even if the collision system under scrutiny 
has not reached high enough energy densities to melt the directly produced J /if) state. 
In particular, it could very well be that the J/ip suppression pattern measured at the 
SPS and RHIC is essentially due to the melting of the if)' and Xc states [3]. 

The picture is made more complex by the fact that already in proton-nucleus 
collisions the charmonium production cross sections scale less than linearly with the 
number of binary nucleon-nucleon collisions. As we will recall in this paper, this 
"normal nuclear absorption" has been seen (at the CERN-SPS and at Fermilab) to 
be significantly stronger for if)' mesons than for J /if) mesons. While the existing Xc 
data are much less accurate, there is no reason to assume that the Xc and J /if) mesons 
have the same "nuclear dependence". A stronger Xc "normal nuclear absorption" 
would decrease the yield of J /if} mesons produced from \c decays and, hence, would 
account for part of the "anomalous J /if) suppression" seen in heavy- ion collisions. 
How much of that "anomaly" might be due to the normal nuclear absorption of the 
ip' and Xc mesons depends on the fractions of J /if) mesons produced by if)' and Xc 
decays. These considerations underline the importance of knowing these fractions, in 
elementary collisions, as accurately as possible. 

In the existing literature, the feed-down fractions are generally assumed to be 
around 10 % for the if)' and around 30 or 40 % for the Xc, usually without mentioning 
experimental measurements or their uncertainties. Yet, the J /if) feed-down fraction 
from if)' decays can be rather precisely determined, from data collected by SPS and 
Fermilab experiments. The Xc case has been much less investigated but recent mea- 
surements, by the HERA-B experiment, indicate a J /if) feed-down fraction from Xc 
decays of around 20 %, considerably lower than the previously assumed values. 

This paper presents a quantitative analysis of the presently available data on feed- 
down contributions to J /if) hadro- production, at fixed target energies. The relevant 
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measurements are presented and reviewed in Section 2. New "world averages" (in- 
cluding uncertainties) of the J/ip feed-down fractions from ip' and Xc decays are then 
derived and discussed, in Sections 3 and 4. 



2 Overview of available measurements 

In this section we briefly review the existing measurements of ip' and \c hadro- 
production, which can be used to constrain the corresponding fractions of indirectly 
produced J/^'s. These fractions are defined with respect to the total (inclusive) J/ip 
yield: 

_ N(J/j> from f) _ a(f) ■ B(f -+ J/j>X) 
[V) iV incl (J/^) ^) ' {) 

and analogous for R(xc)- 

The ip' mainly decays into a J/ip and a pair of pions. However, most experiments 
measure the J/tp and ip' dilepton decays, reporting results for the yield ratio 

*w) ■ bw -» in-) 

This quantity is directly related to the ip'-to-J /ip feed-down fraction, R(ip'), through 
a simple combination of branching ratios, 



pty') = (4.53 ± 0.13) P W) , (3) 



where the numerical values were derived from the PDG tables [4]. 

The R(Xc) values are obtained dividing the number of J/^'s resulting from the 
radiative decays Xc J/4> 7 by the total number of observed J/^'s. 

The running conditions of the experiments providing these measurements are 
summarized in Tables 1 and 2. Most experiments made use of proton or pion beams of 
different energies incident on several target nuclei, but collider experiments have also 
provided some results. The different detectors covered xp intervals extending from 
slightly backward to very forward values. Average xf values have been estimated for 
each experiment, either from the measured distributions or from the variation of the 
statistical errors in the efficiency-corrected spectra. 

While all the relevant experiments are listed in these tables, only a restricted sub- 
sample of data is used in the present analysis. Most importantly, we do not use results 
obtained on the basis of forward xf data. In fact, it is well established (in particular 
by E866 [15], for xp > 0.2) that the J/ip and ip' production cross sections measured 
at high xf, in nuclear targets, exhibit a much stronger nuclear absorption than the 
corresponding mid-rapidity values. A significant role in this behaviour should be 
played by nuclear effects on the parton distribution functions of the target nucleons 
and by the energy loss of the beam partons (or of the produced state) traversing 
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Phase space 


(x F ) 


E331 [51 


D-C 


225 


< Xt? < 7 


~ 0.3 


E444 [61 


D-C 
P 


225 


< ttt < 9 


~ 35 

V7 . - J '.J 


E705 [7] 


p-Li 


300 


-0.1 < x F < 0.5 


~ 0.2 


E288 [81 


n-Rp 


400 


—0 6 <r tu <r 8 


~ 1 


NA38 [9] 


p-W/U 
n-C/Al/Cn/W 

[J V / ill / Kj Li / V V 


200 
450 


-0.4 < y cm < 0.6 


~ 


NA51 rioi 


n-H/R 
p n/ u 


450 


—0 A < ii < 6 


~ 


NA50 96/98 [11] 


n-Rp/Al /On / 

■"■6/ vv 


450 


-0.5 < y cm < 0.5 


~ 


NA50 2000 [12] 


p-Be/Al/Cu/ 
Ag/W/Pb 


400 


-0.425 < y cm < 0.575 


~ 


E771 [13] 


p-Si 


800 


-0.05 <x F < 0.25 


~ 0.1 


E789 [14] 


p-Au 


800 


-0.03 <x F < 0.15 


~ 0.06 


E866 [15] 


p-Re/Fe/W 


800 


-0.1 < x F < 0.8 


~ 0.3 


HERA-B [16] 

L J 


p-C/Ti/W 


920 


-0.35 <x F < 0.1 


-0.065 


WA39 [171 


7T ± -W 


39.5 


-0 5 < x v < 8 


~ 0.2 


E537 [181 


7T~-W 


125 


< Xtt < 1 


~ 0.3 


WA11 [19] 


7r~-Re 


150 


-0.4 <x F < 0.9 


~ 0.3 


E331 [5] 


7T + -C 


225 


< x F < 0.9 


~ 0.35 


E444 [6] 


TT^-C 


225 


< x F < 1 


~ 0.4 


E615 [20] 


7I--W 


253 


0.3 < x F < 1 


~ 0.6 


E705 [7] 


7T ± -Li 


300 


-0.1 < x F < 0.5 


~ 0.2 


E672-706 [21] 


7r _ -Re 


515 


0.1 < x F < 0.8 


~ 0.4 


Experiment 


Collision system 


>/I [GcV] 


Phase space 


(x F ) 


ISR [22] 


PP 


58 (avg.) 


Vcm ^ 






Table 1: Global features characterizing the existing measurements of the ip'-to-J /ip 
cross-section ratio in proton-nucleus, pion-nucleus and proton-proton collisions. 



the nuclear matter. Other effects may also contribute, such as intrinsic heavy-quark 
components of the scattering nucleons and interactions with other produced hadrons 
( "comovers" ) , as discussed in Ref. [33]. The data at high x F may, therefore, reflect a 
non-trivial cocktail of production and absorption mechanisms, certainly not easy to 
disentangle and quantify. For this reason, in this paper we concentrate on the analysis 
of the mid-rapidity data, a choice consistent with our goal of determining reference 
values for the interpretation of the existing observations of quarkonium suppression 
in nucleus-nucleus collisions, also made at mid-rapidity. This means, in particular, 
that the data obtained in pion-nucleus collisions, always significantly extending to- 
wards high x F values, are left out from the present analysis. Further details on the 
data selection are discussed in the following paragraphs. A broader analysis of the 
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Table 2: Global features characterizing the existing measurements of the R(xc) feed- 
down ratio in proton-nucleus, pion-nucleus and proton- (anti)proton collisions. 



existing measurements, trying to take into account possible kinematic dependencies 
induced by nuclear effects, will be the topic of a future investigation. 

The current experimental knowledge concerning ip' production in proton-nucleus 
collisions is essentially determined by the accurate measurements performed by 
NA50/NA51 at the CERN-SPS and by E866 at Fermilab, using several target nuclei. 
Thanks to their several million reconstructed J/ifj events, these measurements provide 
precise determinations of the ip'-to-J/ip cross-section ratio and of its nuclear depen- 
dence. Given their much lower level of precision, the other existing measurements 
(obtained in several kinematical windows and using a single target nucleus) have no 
influence on a global average and were, therefore, left out of the present analysis. 

The NA50 data sets were obtained at 450 [11] and 400 GeV [12], respectively 
with five and six targets. Two statistically independent 450 GeV data samples were 
collected, at different proton beam intensities. In our study we used their average, 
after considering the corrections reported in Ref. [12]. 

The J /ip and ip' results of E866 [15] were reported as ratios between the yields 
obtained with heavy and light targets (W/Be and Fe/Be), as a function of x F . These 
measurements provide heavy-over-light ratios of p(i/)'), which cannot help determin- 
ing the feed-down fraction value but constrain the difference between the nuclear 
absorption rates of the two charmonium states. Although data points exist up to 
x F = 0.8, we restricted our study to the range x F < 0.2, where only the W/Be ratio 
was reported. 

Also the Xc measurements mentioned in Table 2 deserve a few remarks. The error 
we quote for the E705 value reflects the systematic uncertainties mentioned in their 
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paper [24]. Although the E771 publication [25] provides no explicit value for R(x c ), 
it should be possible to derive it from the quoted Xci and Xc2 cross sections or, alter- 
natively, from the yields of reconstructed Xci, Xc2 and J/ip mesons, using the quoted 
efficiencies. It turns out, however, that these two methods lead to significantly dif- 
ferent results. Moreover, the information provided is insufficient to properly evaluate 
the R{Xc) uncertainty. Therefore, we did not consider the E771 measurement in our 
analysis. 

It is worth noting that the HERA-B 2003 R(Xc) results [27] include a systematic 
uncertainty (of around 10 %) due to the dependence of the detector's acceptance on 
the assumed J/ip polarization, taking into account that J/^'s from Xc decays may 
have a polarisation different from the directly produced ones. This effect was not 
considered by the previous experiments, given the poor statistical accuracy of their 
measurements. 

Given the considerations expressed above, we have selected for our analysis the 
pijp') and R(xc) measurements listed in Table 3. 

3 J/i/j feed-down from -0' decays 

The experimental points selected for the determination of the J/ijj feed-down con- 
tribution from i])' decays are shown in Fig. 1, as a function of the size of the target 
nucleus (in the case of the SPS data) or of xp (in the case of E866). These measure- 
ments clearly show that the ip' and J/tp states are differently absorbed by the nuclear 
medium. 

In order to determine the ip'-to-J /i/j feed-down fraction in pp collisions, R (ip'), 
all selected measurements were simultaneously fitted within the framework of the 
Glauber formalism [34] , using the so-called u pL parametrization" : 

cr(pA -> ip) I Aa(pN -> i/j) = exp(-o- abs pL) , 

where p is the nuclear density and L is the nuclear path length traversed by the 
charmonium state, of absorption cross section <j a ^ s . The pL values were determined 
through a Glauber calculation, for each nuclear target, taking into account the ap- 
propriate nuclear density profiles, as described in Ref. [12]. The fit provides two 
parameters: the R°(ip') "reference" feed-down fraction (corresponding to L — 0) and 
the difference between the ip' and J/ip absorption cross sections, where the J /if) term 
does not include the ijj' decay contribution (to remove auto-correlation effects). 

It should be kept in mind that this parametrization represents a rather simplified 
description of the nuclear absorption process, convoluting in a single effective absorp- 
tion cross section a multitude of physical effects. In particular, <7 a b s is assumed to 
be a "universal quantity" , independent of the collision energy and of the kinematical 
properties of the produced charmonium states. It also implicitly incorporates the 
nuclear modifications of the parton distribution functions, possible energy loss mech- 
anisms, formation time effects, etc. Following most previous studies of charmonium 
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Table 3: The p(ip') and R(xc) measurements selected for the present analysis. The 
L values correspond to an average nuclear density of 0.17 fnT~ 3 . 



absorption in nuclear matter, we use this simple parametrization in the analysis pre- 
sented in this paper, where we focus on the mid-rapidity results; this issue will be 
revisited in the future, in the scope of a broader investigation. 

A global fit to all data points leads to the dashed lines in Fig. 1, with a chi-square 
probability of only 1 %, clearly indicating that the model is unable to properly account 
for the NA51 measurements, performed with hydrogen and deuterium targets (the 
two leftmost points in the middle panel). Maybe the fact that protons and deuterons 
are exceptionally light nuclei places them out of the domain of applicability of the 
model we are using because they are not large enough to be traversed by fully formed 
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Figure 1: The p(ip') values as a function of the nuclear path length, L, from NA50/51 
measurements at 400 GeV (a) and 450 GeV (b), and the p(^')w / p(ip')Be ratio mea- 
sured by E866, in the —0.1 < xp < 0.2 window (c). The curves are the result of the 
global fit described in the text, including (dashed lines) or excluding (solid lines) the 
NA51 points. 



charmonium states. It should also be noted that the use of "nuclear density profiles" 
in the Glauber calculation of the proton and deuteron p L values is not as reliable as 
in the case of the heavier nuclei. Furthermore, it is not clear that the same value, 
0.17 nucleon/fm 3 , should be used as average nuclear density for all nuclei, including 
protons and deuterons, when extracting L from the calculated pL values. Without 
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Figure 2: 68 % (solid lines) and 99 % (dashed lines) confidence level contours for 
the bi-dimensional probability distribution of the parameters R°(ip') and cr a b s (V' / ) — 
cr abs(J/V ; + Xc)- The thick black contours delimit the region favoured by the global 
fit, while the thin coloured (or grey) ones reflect the individual data sets. 



the pp and p-D points, the best description of the data is represented by the solid lines, 
with a chi-square probability of 27 %, reflecting a much better compatibility between 
the data and the model used in our fit. The corresponding feed-down fraction is 

i?V) = (8-1 ±0.3)% . (4) 

Including the NA51 points decreases the result to (7.9 ± 0.3) %, a negligible change 
despite the visible degradation of the fit quality. 

The correlation between the two fit parameters is shown in Fig. 2, as a bi- 
dimensional contour plot. Although the three data sets give compatible results, they 
nevertheless indicate that the difference between the charmonium absorption cross 
sections decreases with increasing collision energy. 



4 J/V* feed-down from % c decays 

The R{Xc) values collected in Table 3 are shown in Fig. 3 as a function of L. The 
curve is the result of a fit analogous to the one explained in the previous section, 
using the "pL parametrization" and leaving free the difference between the effective 
absorption cross sections of the \ c mesons and of the J/ip mesons not coming from 
Xc decays. Given the conjecture, suggested by the if)' analysis, that measurements 
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Figure 3: The R(Xc) measurements used in the present analysis as a function of 
the nuclear path length L. The curve is the result of the fit described in the text 
(excluding the first point). 



performed with very light nuclei are not accountable within the simple absorption 
model adopted here, the pp point is excluded from the fit. The resulting feed-down 
fraction (for L = 0) is 

R\ X c) = (25 ± 5) % , (5) 

with a fit x 2 probability of 25 %. 

The R°(xc) value considerably depends on the difference between the absorption 
cross sections of the two charmonium states (see Fig. 4). Therefore, a more precise 
R°(Xc) value can be obtained if an improved understanding of charmonium absorption 
in nuclear targets significantly reduces the allowed range of <J a b s (Xc) — o"abs(J/V')- 



5 Summary 

We presented and reviewed the presently available ip' and Xc hadro-pro duct ion mea- 
surements, and derived global averages of the J/ip feed-down fractions from ip' and 
Xc decays, at mid-rapidity: 

#V) = (8-1 ± 0-3) % , R°(Xc) = (25 ±5)% . (6) 

These averages reflect measurements performed at collision energies up to y/s ~ 
60 GeV. At much higher energies, CDF measured R(xc) = (30 ± 7) % in pp collisions 
at y/s = 1800 GeV [32] and PHENIX reported preliminary values obtained in pp 
collisions at y/S = 200 GeV: R(ip') = (8.6 ± 2.5)% and R(x c ) < 42% (at 90% 
C.L.) [35]. More precise measurements would be needed to probe an eventual energy 
dependence of the J ftp feed-down fractions from decays of heavier charmonium states. 

Since most of the existing measurements were performed with nuclear targets, the 
derivation of the R°(i/j ; ) and R°(xc) values relevant for elementary collisions (L = 0) 
requires modelling the influence of the nuclei on the production yields of the three 
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Figure 4: 68 % and 99 % confidence level contours for the bi-dimensional probability 
distribution of the fit parameters R°(Xc) and d a b s (Xc) — Oabsp/^ + ip')- 

charmonium states. In the study reported in this paper we followed a widespread 
model where the three charmonium states are analogously absorbed while traversing 
the nuclear matter, with three (a priori different) effective absorption rates. This 
rather simple model provides a reasonable description of the available measurements 
if we restrict the analysis to the mid-rapidity data and exclude values obtained with 
"exceptionally light nuclei" (protons and deuterons) . The extension of our analysis to 
a broader set of measurements requires an improved phenomenological model, which 
should reflect the following observations. It is only for nuclear targets heavier than 
beryllium that the J/tp and ip' nuclear absorption rates are significantly different. 
This difference decreases when the collision energy increases and when we approach 
forward xf. The R{Xc) value derived from data collected with heavy nuclei and 
at small |xf|, 0.22 ± 0.03, is significantly smaller than the value derived from data 
collected with light nuclei and at forward xp, 0.36 ± 0.02. 

These might be indications that a proper understanding of charmonium absorp- 
tion requires considering that the objects traversing the nuclear matter are not fully 
formed J/ip, ip' or Xc states but rather pre-resonance states having a suitable time 
evolution. This is the topic of a more complex investigation, to be reported in a 
future publication. 
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work was partially supported by the Fundagao para a Ciencia e a Tecnologia, 
Portugal, under contracts SFRH/BPD/42343/2007, SFRH/BPD/42138/2007 and 
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